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Abstract. In this paper we introduce kinetic equations for the evolution of 
the probability distribution of two goods among a huge population of agents. 
The leading idea is to describe the trading of these goods by means of some 
fundamental rules in price theory, in particular by using Cobb-Douglas utility 
functions for the binary exchange, and the Edgeworth box for the description 
of the common exchange area in which utility is increasing for both agents. 
This leads to a Boltzmann-type equation in which the post-interaction vari- 
ables depend in a nonlinear way from the pre-interaction ones. Other models 
will be derived, by suitably linearizing this Boltzmann equation. In presence 
of uncertainty in the exchanges, it is shown that the solution to some of the 
linearized kinetic equations develop Pareto tails, where the Pareto index de- 
pends on the ratio between the gain and the variance of the uncertainty. In 
particular, the result holds true for the solution of a drift-diffusion equation 
of Fokker-Planck type, obtained from the linear Boltzmann equation as the 
limit of quasi-invariant trades. 

Keywords. Wealth and income distributions, kinetic models, Boltzmann 
equation, Fokker-Planck equation. 



In an effort to understand the emergent properties appearing in complex agent- 
based systems, various concepts and techniques of statistical mechanics have been 
fruitfully applied in the last twenty years to both social and economic fields. This 
is well documented both by recent books on these topics [3J [T^l HH1 S3 113 US] , as 
well as by the introductory articles [THl HH HH [301 HH]- in particular, the point of 
view of collisional kinetic theory revealed to be a powerful instrument to describe 
the behavior of systems in which the mechanism of variation could be described 
mostly in terms of binary interactions between agents [17l [TBI El Ell EI] ■ In kinetic 
theory of rarefied gas, where binary collisions between molecules are the dominant 
phenomenon, the evolution of the (spatially uniform) density of the molecules is 
described by the spatially homogeneous Boltzmann equation [T^l [T2j , 



In the bilinear Q term accounts for all kinematically possible (those that con- 
serve both momentum and energy) binary collisions. The post-collisional velocities 
(v*,w*) are consequently given by the linear (in terms of the pre-collisional ones 
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1. Introduction 
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(u, w)) relations 

(1.2) v* — —(v + w + \v — w|n), w* = —(v + w — \v — w\n), 

in which n is a unit vector. The Boltzmann equation (jl.ip describes the relaxation 
of the gas towards equilibrium, which in this case is given by the Maxwellian density 
[l2l H3] of mass p, velocity U and temperature T, that is 

(1.3) M{v):= P^ exp {Jl^l\. 
y ' W (2ttT) 3 / 2 P \ 2T J 

One of the most used models for the collisional operator Q is the so-called Maxwellian 
molecules operator [2j EJ El [8] , given by 

(1.4) Q(f)(v)= [ B((v-w)-n)[f(v*)g(w*)}dwdn-pf(v). 

JR 3 xS 2 

In p. 41) £?(•) is a measure of the collision frequency, dn denotes the normalized 
surface measure on the unit sphere S 2 , and 

P= f(v) dv . 

For Maxwellian molecules the operator simplifies, and the positive part of it takes 
the form of a generalized convolution [7]. 

Boltzmann-likc equations for the evolution in time of the density of wealth in 
a market composed by many agents are easily recovered by identifying molecules 
with agents, particles 'energy with agents' wealth, and binary collisions with trade 
interactions. In this way, the kinetic description of market models via a Boltzmann- 
type equation provides one possible explanation for the development of universal 
profiles in wealth distributions of real economies. 

To obtain the kinetic description one usually resorts to few reasonable assump- 
tions, that date back to the work of Angle pQ. First, agents are considered indis- 
tinguishable. Then, an agent's state at any instant of time t > is completely 
characterized by his current wealth w > 0. When two agents encounter in a trade, 
their pre-trade wealths (v,w) change into the post-trade wealths (v*,w*) according 
to the (linear) rule 

(1.5) v* = p\V + qiw, w* — q2V + P2W . 

In most of the existing models, the interaction coefficients pi and qi are non-negative 
random variables. While q± denotes the fraction of the second agent's wealth trans- 
ferred to the first agent, the difference p\ — q^ is the relative gain (or loss) of wealth 
of the first agent due to market risks. One generally assumes that pi and q{ have 
fixed laws, which are independent of v and w, and of time. 

In these one-dimensional models, the probability distribution f(w,t) of wealth 
of the ensemble coincides with agent density and satisfies the associated spatially 
homogeneous Boltzmann equation of Maxwell type (|1.4|) 



(1.6) 



-£ = Q+(f,f)-f, 



KINETIC MODELS FOR THE TRADING OF GOODS 



3 



on the real half line, w > 0. The collisional gain operator Q+(f, /), which quantifies 
the gain of wealth v at time t due to binary trades, acts on test functions tp(w) as 



with (•) denoting the expectation with respect to the random coefficients pi and 
qi in (|1.5p . The large-time behavior of the density is heavily dependent of the 
evolution of the average wealth [24l [39] 



Conservative models are such that the average wealth of the society is conserved 
with time, M(t) — M, where the value of M is finite. In terms of the interaction 
coefficients, this is equivalent to impose the conditions (p± + (72) = (P2 + <Zi) = 1 • 

The Boltzmann equation p. 61) belongs to the Maxwell type. As in operator fll.40 
the collision frequency is independent of the relative velocity, and the loss term in 
the collision operator is linear. This introduces a great simplification, that allows 
to use most of the well established techniques developed for the three-dimensional 
spatially homogeneous Boltzmann equation for Maxwell molecules in the field of 
wealth redistribution [2TJ HI HQ 0S2 [Ml \M \M ■ 

In general, the richness of the steady states for kinetic market models is the main 
remarkable difference to the theory of Maxwell molecules [7]. While the Maxwell 
distribution (|1.3j) is the universal steady profile for the velocity distribution of 
molecular gases, the stationary profiles for wealth can be manifold, and are in gen- 
eral not explicitly known analytically. In fact, they depend heavily on the precise 
form of the microscopic modeling of trade interactions. Consequently, in investiga- 
tions of the large-time behavior of the wealth distribution, one is typically limited 
to describe a few analytically accessible properties (e.g. moments and smoothness) 
of the latter. 

This richness of steady states makes both the theoretical and numerical study of 
these models highly interesting. Indeed, it corresponds mathematically to identify 
the limit distribution in a problem which has many analogies with a generalized 
central limit theorem. Also, it allows to clarify the types of binary interactions 
which produce fat tails in the distribution at infinity (Pareto tails), which are 
shown to form in the profiles of the distribution of wealth in western societies. 

On the other hand, the socio-economic behavior of a (real) population of agents 
is extremely complex. Apart from elements from mathematics and economics, a 
sound description — if one at all exists — would necessarily need contributions from 
various other fields, including psychology. Clearly, the majority of the available 
mathematical models are too simple to even pretend to reflect part of the real 
situation. 

For a better understanding of some of the main outcomes of markets in western 
societies, in recent years the kinetic community started to improve the economic 
modeling of agent systems by including more realistic assumptions, without making 
too heavy the underlying Boltzmann equations. Two papers in particular fit into 
this line of thinking. 



(1.7) 




(1.8) 
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In [20] Cordier, Pareschi and Piatecki derive a kinetic description of the behav- 
ior of a simple financial market where each agent can create their own portfolio 
between two investment alternatives: a stock and a bond. In this case, the varia- 
tion of density is not based on linear binary collisions like (jl.5[) . but it is derived 
starting from the more realistic Levy-Levy-Solomon microscopic model for price 
formation [HH [35] . The model in [20] consists of a linear Boltzmann equation for 
the wealth distribution of the agents coupled with an equation for the price of the 
stock. From this model, under a suitable scaling, Cordier, Pareschi and Piatecki 
derive a Fokker-Planck equation and show that the equation admits a self-similar 
lognormal behavior. For the first time, the kinetic model in [20j attempts to join to 
simple financial rules a kinetic equation of Boltzmann type, able to describe a com- 
plex behavior that could then mimic the market and explain the price formation 
mechanism. 

A second interesting example of this coupling has been recently proposed in 
[38] by Maldarella and Pareschi. They introduce a relatively simple model for a 
financial market characterized by a single stock or good and an interplay between 
two different trader populations, chartists and fundamentalists, which determine 
the price dynamics of the stock. The model has been inspired by the microscopic 
Lux-Marchesi model |36[ I37j . The main novelty here is to couple the financial rules 
with the opinion of traders through a kinetic model of opinion formation recently 
introduced in [47] . Moreover, some psychological and behavioral components of 
the agents, like the way they interact with each other and perceive the risk, which 
may produce nonrational behaviors, are taken into account. This is done by means 
of a suitable value function in agreement with the Prospect Theory by Kahncman 
and Tversky [321 133] . As they show, people systematically overreacting produces 
substantial instabilities in the stock market. 

These two papers indicate that the methods of kinetic theory can be fruitfully 
used to understand the market behavior even in presence of more realistic ways 
of interactions between agents. The objective weakness of models based on linear 
binary trades of type (jl.5[) is linked to the fact that people trades according to rules 
which are in general difficult to justify from a microeconomic point of view. The 
common ingredients which appear in these trades are the saving propensity concept 
[14] (agents usually do not trade with all their wealth), and the risk concept [21] 
(there is no certainty of gain). While highly reasonable, this type of interaction 
does not clarify while agents are pushed to trade. 

One of the fundamental assumption of prize theory in economics is that people 
trades to improve its utility. Among binary interactions which are motivated by an 
increasing of the individual utility, one of the most used pictures is furnished by the 
Edgeworth box [26] , which is frequently used in general equilibrium theory, and can 
aid in representing the competitive equilibrium of a simple system or a range of such 
outcomes that satisfy economic efficiency. Edgeworth box can fruitfully be applied 
in presence of an agent-based system in which agents possess a finite number of 
goods of n > 2 different types. Inspired by this mechanism of increasing utility 
and competitive equilibrium, we will introduce in the following a kinetic equation 
of Boltzmann type for the evolution of the distribution density of the quantities of 
two goods, said (x,y), in a system of agents. In other words, we will be concerned 
with the time evolution of the probability density f(x,y,t) of having x goods of 
the first type and y goods of the second type at time t > 0, given a known initial 
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distribution of these goods at time t = 0. The variation in time of f(x,y,t) will 
be entirely due to the microscopic binary interaction driven by the Edgeworth box 
structure. We will take into account the incomplete knowledge of the game in the 
binary interaction by allowing agents to be wrong with respect to the ideal outcome 
predicted by the increase of utility. This will be done by introducing a randomness 
in the Edgeworth box outcome. 

In more details, in the next Section we will resume the Edgeworth box interac- 
tion, and the consequent Boltzmann-type equation, together with its main prop- 
erties. Surprisingly, on the contrary to most of the previous types of trade, this 
exchange rule leads to a highly nonlinear binary interaction, which is difficult to 
handle, if not numerically. For this reason, in Section [3] we will resort to a suit- 
able linear Boltzmann equation, which is obtained by allowing the agent to interact 
(according to Edgeworth box), simultaneously with a sufficiently high number of 
agents in the market. By Fourier based methods, it will be shown that, under 
reasonable conditions on the initial data, this linear equation has a unique solution, 
and the steady states are concentrated along a well-defined line. Finally, in Section 
IU we will resort to an asymptotic procedure (the so-called quasi invariant trade 
limit), to obtain a linear Fokker-Planck equation which describes the essentials of 
the trading of goods, and it is relatively easy to treat from a mathematical point of 
view. In particular, it will be shown that the solution converges towards a steady 
state with fat tails. 



As discussed in the Introduction, most of the existing kinetic models for wealth 
distribution are based on rigid assumptions which, if on one hand can be shared, 
from the other hand are not deeply related to economic principles, like price theory 
[27j . The aim of this Section is to introduce a new framework for trades, which is 
derived directly from the basic principles of economy. 

Individuals exchange goods. The benefits they receive depend on how much 
they exchange and on what terms. Price theory tries to give an answer to this 
fundamental question. In the case of a binary trade, there may be many different 
exchanges, each of which would be beneficial to both parties; some exchanges will be 
preferred by one person, some by the other. There are then two different questions 
to be settled. One is how to squeeze as much total gain as possible out of the 
opportunities for trade; the other is how that gain is to be divided. The two 
individuals who are trading have a common interest in getting as much total gain 
as possible but are likely to disagree about the division. Let us suppose in the 
following that agents in the system possess the same two types of goods, and there 
is no production, so that the total amount of goods remains unchanged. If (xa, Ha) 
(( x b,Vb) respectively) denote the quantity of goods of two agents A and B, the 
quantities 



2. The Edgeworth box 



(2.9) 



VA 



PA = 



XA + Xb 1 



qA = 



VA + VB 



are the percentages of goods of the first agent. Note that the point (pa, qA) belongs 
to the square S — [0, 1] x [0, 1]. In order to give a meaning to the reasons of trading, 
it is classical to assume that agent's behavior is driven by a utility function. One 
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of the most popular of these functions is the Cobb-Douglas utility function 
(2.10) U( P ,q)=p a qP, a + 0=1. 

Each agent will tend to maximize its utility by trading. The values a and /3 are 
linked to the preferences that the agent assigns to the two goods. If a > (3, the agent 
prefers to possess goods of the first type (numbered by x). The choice a = f3 = 1/2 
clearly means that the two goods are equally important for A. Given the percentage 
point (pa, 5a) of the agent A, the curve 

U(p,q) = U(p A ,qA) 

denotes the indifference curve for the agent A. Indeed, any point on the indifference 
curve has for A the same utility. Note that the indifference curve for A entirely 
belongs to <S, and splits the square into the two regions 

U A = {{p,q) : U(p,q) < U{p A ,q A )}, U\ = {(p,q) : U{p,q) > U{p A ,q A )}. 

Clearly, any trade which will move the percentages of agent A into the region U\ 
will increase its utility function, and will be acceptable for the agent. 

In the case of a binary trade, also agent B has a Cobb-Douglas utility function, 
with parameters of preference which in general are different from the parameters of 
agent A. Likewise, it will be an indifference curve for B, and a good region in which 
the utility of B is increased after trading. Ultimately, a trade will be acceptable 
for both agents, if their percentages after the trade belong to the regions in which 
their utilities are increased. 

An ingenious way of looking at such a situation is the Edgeworth Box, named 
after Francis Y. Edgeworth, the author of a nineteenth century work on economics 
called Mathematical Psychics [26] . 

The Edgeworth box is built up by rotating the square S in which is designed the 
indifference curve of agent B of 180° around the center of the square (1/2, 1/2), 
and considering together the indifference curve of A and the rotated indifference 
curve of B on the same square. 

Any point inside the box represents a possible division of the total quantity of 
A and B, with the share of A measured from the lower left-hand corner, and the 
share of B from the upper right-hand corner. Any possible trade is represented by 
a movement from one point in the box, to another. 

By construction, since the Cobb-Douglas utility function is jointly convex, the 
utility of A increases if the point moves up and to the right; the utility of B increases 
if the point moves down and to the left. The region in which this happens is exactly 
the region common to U\ and to Ug rotated of 180°. 

This operative way of trading furnished by the Edgeworth box, has interesting 
consequences. If the percentages of A and B are such that there is a possible region 
for trading, once the trade has been done, the point has moved inside the region, 
and the new region for trading is smaller. If the two agents were so smart to choose 
a point inside the region in which the two indifference curves are tangent each other, 
since they curve in opposite directions, this means that starting from this point, 
any point that is on a higher indifference curve for the first agent must be on a 
lower curve for the second; any trade that makes one better off makes the other 
worse off. The point in which the two curves are tangent each other is not unique. 
The set of all points from which no further mutually beneficial trading is possible 
(Pareto optimal points) is called the contract curve. 
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The contract curve is related to prices. At that tangency of the two indifference 
curves, in fact, the slope of the tangency line represents the relative prices for the 
two goods. Hence, there are relative prices that will be consistent with the Pareto 
optimum, and these prices will maximize the possible budget for both agents. 

This short discussion on utility functions, trading and prices is the basis of our 
construction of the binary trade. In what follows, we suppose that all agents in the 
system have the same utility function (|2.10|) . If agent A has percentages p and q 
in his Hedgeworth box, we consider as possible trades the movements of the point 
(p, q) into (p*,q*) € S, where 

p* =p + \j3(q-p) +n[q-p) 

q* =q+\a(p-q) + fl(p-q). 

In (12. 11|1 < A < 1, and /i and p, are random variables with zero mean and finite 
variance. We will assume moreover that 

(2.12) 0<A/3 + /x<l, 0<\a + p, < 1. 

Under this condition, the post-trade point (p*,q*) belongs to S, and it is an ad- 
missible point for the Edgeworth box. It is immediate to show that, unless p = q, 
and in absence of randomness, that is if 

(2 13) p* = p + \!3(q - p) 

q* = q + \a(p-q), 
the trade increases the Cobb-Douglas utility function of agent A. Indeed 

^■U( P *,q*) = a/J(l - A)(p - g)V-V _1 > 0, 
which, coupled with the convexity of the utility function gives 

U(p*,q*) = Ufa q) + a/3( P - g) V" V _1 A + \^ U $* ^ 

U{p, q) + ap{p - g) 2 p a_ V _1 A. 
Note moreover that the trade (|2.13l) is such that p* ^ q* , unless p — q. Concerning 
agent B, since it starts from the point (1 — p, 1 — q) E S, according to (j2 . 13[) it 
moves to the point 

p=l-p + A0(l-g-(l-p)) 
(2.14) _ 

q=l-q + Aa(l - p - (1 - g)), 

which, proceeding as before, implies 

U(p,q) > U(l-p,l-q). 

Hence, in absence of randomness, the trade (|2.13l) is such that the Cobb-Douglas 
utility functions of both agents are increased, locating the new point in S in the 
interior of the allowed area for trades predicted by the Edgeworth box. Clearly, 
this trade refers to the purely theoretical situation in which the traders A and B 
know all about the trade and its result. In general, this is not realistic, and what 
traders can hope, after a certain number of trades, is that the majority of their 
trades had the result of increasing their utility. This can be expressed by saying 
that a reasonable result for traders is to obtain, in a single trade, the increasing of 
their utility only in the mean. If one agrees with this, then trades of type (|2. 1 1|) 
satisfy all constraints of Edgeworth box, but the single trade can move the point 
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in a region which is convenient only for one of the two traders, as well as in the 
region which is not convenient for both. 

Once the trade has been done, according to the rules of price theory, the quan- 
tities of goods of agents A and B have been changed accordingly. Making use of 
the constraints xa + %b — x* A + x* B and ua + Ub — Ua^~ Ub^ P-HD implies for the 
agents A and B the post-trade quantities of goods 



We remark that, on the contrary to what is usually assumed in the wealth trades of 
type (|1.5[) . the post-trade quantities in (|2.15[) are related to the pre-trade quantities 
by nonlinear relations. It can be verified, however, that the constraints about the 
conservation of the total quantities of goods in the trade are verified, and 

x* A + x B =x A + x B , y* A + y* B = ija + Vb- 
3. A LINEAR BOLTZMANN EQUATION FOR TRADING OF GOODS 

Let f(x,y,t) denote the density of agents with quantities x and y of the two 
goods at time t > 0. Without loss of generality, we will assume in the following 
that (x, y) are nonnegative real numbers. A Boltzmann-like equation of Maxwell 
type for the time evolution of f(x, y, t) can be written by observing that the collision 
rules are here given by (|2. 15|) . A useful way of writing the Boltzmann equation has 
been recalled in the introduction, where the Boltzmann equation (|1.6j) has been 
rewritten resorting to the weak form (|1.7p . It corresponds to write, for any given 
smooth function ip, the spatially homogeneous Boltzmann equation in the form 



The right-hand side of equation (|3.16[) describes the change of density due to colli- 
sion of type (|2.15[) . Given two agents with quantities of goods given by (x,y) and 
(x 1,2/1), one obtains from (j2.15[) the post-trade quantities 



(2.15) 




(3.16) 




(3.17) 




Hence 



(3.18) 





(<f(x*,y*) - <p{x,y))f{x,y,t)f(x 1 ,y 1 ,t)dxdydx 1 dy 1 
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In alternative, one can use the identity 
(3.19) 



l f{x,y)Q(f,f)dxdy = -( [ {<p(x*,y* + <p(xl,yfl 



2 



-Vfav) ~ l f{xi,yi))f(x,y 1 t)f(xi 1 yi,t)dxdydxidy l 
By choosing ip(x, y) = x (respectively <p(x, y) — y) one verifies that the mean values 

m x (t) = I xf(x,y,t)dxdy 



(3.20) 

™>v(t) = I yf(x,y,t)dxdy, 



remain constant in time m x (t) = m x (0) (respectively m y (t) — m y (0) ). These 
conservations are consequence of the constraints of conservation of the quantities of 
goods in the trades. Since the exchanges of goods of type (|3.17l) are nonlinear, the 
study of the properties of the solution to the Boltzmann equation (|3.16p is difficult. 
In particular, it is cumbersome to find the evolution of the higher moments of the 
solution in a closed form. Hence, the Boltmann equation (|3.16p is the starting 
point for a numerical study of the evolution of the density by means of MonteCarlo 
methods [9j|42]. 

In order to derive a more treatable model, which maintains most of the features 
of the nonlinear one described by equation (|3.16[) . let us study in more details 
the exchange rule driven by the Edgeworth box, assuming that, instead of single 
binary trades, any agent is allowed to exchange goods with part of the agent market, 
composed by a certain number of agents. We remark that the exchange rule (|2.15[) 
for the quantity of the first good can be fruitfully rewritten as 
(3.21) 

I XA + Xb \ ( XAJrXB \ 

x* A = x A + (A/3 + n) y yA + y ^ yA - x aJ =x a + (X/3 + fi) f yA 2 y A - x A j ■ 

This writing clarifies that, in the post-trade quantity of goods in (|3.2ip . the coef- 
ficient of yA is the ratio between the mean values of the quantities x and y of the 
two agents exchanging goods. 

If one agent exchanges goods, according to the Edgeworth box strategy, with a 
number of N agents, the quantity of the first good changes according to 

. a . n / X A + Xi -\ h x N 

x A = x A + + (i) ■ ■ ■ y A - x A 

\VA+yi-\ h Vn 



. and, if N is sufficiently large 
(3.22) 



xa + xi H + x N _ ]v+T ^ m x 



VA + yi-i \-yN 



VA+yi-i hyjv 

JV+1 



Since the ratio between the mean values is a conserved quantity, by using approx- 
imation (I3.22[) we obtain a consistent linear version of the trade of goods, clearly 
valid under the assumptions specified above. Therefore, by allowing agents to in- 
teract with part of the market, according to the rules of the Edgeworth box, the 
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in which the action of the market appears only through the mean values m x and 
rriy. Consequently, if f(x, y, t) denotes the density of agents with quantities x and y 
of the two goods at time t > 0, f(x, y, t) satisfies now the linear Boltzmann equation 



Also for equation (|3.24j) . it can be easily verified that the mean values (|3.20|) remain 
constant in time. 

Due to its linearity, the Boltzmann equation (|3.24p can be studied by resorting to 
Fourier transform, by standard methods available for the Boltzmann equation for 
Maxwell molecules, and their extensions to kinetic models of wealth distribution [24j 
[39] . Among other properties, the analysis of [39] in terms of Fourier-based metrics, 
allows to prove existence (and uniqueness) of a steady distribution f oa (x, y), under 
few reasonable conditions on the initial distribution of goods, like existence of a 
certain number of principal moments. 

A quite useful expression for the Fourier transformed equation follows by writing 
the linear interaction (|3.23[) as 

m,yX* = rriyX + (A/3 + /i) (m x y — m y x) 

m x y — m x y + (Xa + fi) [m y X - m x y) . 

The trade (|3. 25|) corresponds to a dissipative interaction of particles with a random 
addition proportional to the relative velocity, similar to the one studied in |10) . 
Then, the pair (x' = m y x, y' = m x y) can be fruitfully used as new set of variables, 
and f(x',y',t) still satisfies the Boltzmann equation (|3.24p . By choosing 

ip(x, y) = exp {-i(x£ + yrf)} , 

and using (|3.25j) . one directly obtain from (|3.24[) the equation for the Fourier trans- 
form / of /, that takes the form 



Note that, by conditions (12. f 2j) . the random coefficients in (I3.27|) are non negative. 
It is immediate to check that the post-trade Fourier variables satisfy the constraint 



(3.24) 




(3-26) -1MU. -fc v ). 
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< 



Hence, the solution to the Boltzmann equation (|3.24p will concentrate for large 
times along the line m x y — m y x. 

To verify that this behavior holds true, one resorts to previous works on similar 
subjects. The standard Fourier-based norm to be used in connection with this type 
of equations [24] is 

3.30 d s (f,g) = sup , , s > 0. 

s, V jto (\£,r + \vn s/2 

Relationships of this norm with other equivalent norms widely used in probability 
theory and mass transportation can be found in 11 . Uniqueness follows by the 
following argument. Let fi(£,T],t) and fi(£,T],t) two solutions to equation (|3.26[) 
corresponding to the initial data /o,i(£,?7) and fo,2(£,v)- Then it holds 

d M,y,t)-f 2 (t,T,,t) M,y,t)-f 2 ^,y,t) _ 

at (iei 2 + |?7l 2 ) s/2 (l£l 2 + M 2 ) s/2 
/ fi(C,y*,t) - M*,v*,t) (iri 2 + w\ 2 y^ 

\ (iei 2 + hi 2 ) s/2 m 2 + \v\ 2 ) s/2 , 
d( f f]sm / (\e\ 2 + \vrr /2 * 

Suppose that 

,,,,, r / (in 2 + hi 2 ) s / 2 \ , , 

(3 - 31) Cs= ^\lW/ <+M 

Then, application of Gronwall's inequality gives (cfr. |28j ) 

(3.32) d s (h(t)J 2 (t)) < 4(/o,i,/o,2)exp{(C s - 1))*}. 

From (|3.32p , uniqueness of solution in a fixed time interval follows. The case in 
which C s < 1 leads to better decay. In fact, if this is the case, for all times t > the 
d s -distance of f± and f 2 decays exponentially in time, and, proceeding as in [39] it 
can be proven that, if f(x, y, t) is a weak solution of the Boltzmann equation (|3.24l) . 
which has initially finite moments up to order 2, then / converges exponentially 
fast in d s to a steady state foo . In addition has the same mean of the initial 
datum, and it is the only steady state with this mean wealth. 

To verify that the constant C' s can be strictly less than one, it is in general not 
direct. However, the following argument shows that we can easily prove that the 
<i s -norm is at least non increasing. Let us set 

A = A/3 + n, B = Xa + fi. 

In reason of condition (|2.12[) . both A and B are strictly positive. It is immediate 
to recognize that 

(3.33) AC -Brf = (1 - A - B)(A$ - Brf) 

Therefore, if the random variables /i and fl, in addition to (|2.12p satisfy the addi- 
tional constraint 

(3.34) (\l-A-B\ s ) = (\l-\-ti-ii\ s ) <l 
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(\C+V*\ 2 + \M*-By*\ 2 r 

(\a + v\ 2 + \M-B V \zy 



) 



< 1, 



Hence, by passing to the new variables 



a = f + v, m = M~ Br] 



we can prove that, with respect to this set of variables the d s (f,g) is not increasing 
for all s > 0. This argument implies the global existence of a unique solution to the 
Boltzmann equation p.24j) . but it is not conclusive with respect to the convergence 
of the global solution towards a unique (in terms of the initial density) steady 
state. Also, it is not clear if, along the line in which the measure steady solution 
is concentrated, there is a behavior at infinity with fat tails. This is related to the 
possibility that, under the constraint of conservation of the total number of goods 
in a binary exchange, there is a possibility to form a class of rich agents. 

Instead of working directly on the Boltzmann equation (|3.24p . in the next Section 
we will deal with the derivation of a Fokker-Planck equation which, while keeping 
the main properties of the Boltzmann equation, will be easier to treat in connection 
with the aforementioned question about tails. 



As often happens when dealing with models of multi-agent systems, one of the 
fundamental issues is to understand the emergent properties appearing in reason 
of the type of interactions. In the problem under study, this corresponds to the 
knowledge of the properties of the equilibrium distribution foa(x,y). To this aim, 
let us introduce a new set of variables, which are more adapted to our purposes. 
Starting from (|3.25p . let us define 

(4.35) v = m y x + m x y, w — m y x — m x y. 

Note that v € R+, while w E R. In addition, by construction < v. With respect 
to the new variables (v, w), the trade induced by the Edgeworth box reads 



These new variables better enlighten the outcome of the exchange, and the con- 
sequences we described by Fourier transform methods in Section [3J In absence of 
randomness, the interaction is dissipative in the w variable, since 



The same property remains true if the random variables /i and jl are such that, for 
a given A 



Note that condition (|4.37[) is the analogous of condition (I3.34[) . which relates the 
smallness of the support of the random variables to the size of A. If (|4.37[) holds 
true, in view of the dissipation, the solution to the linear Boltzmann equation (|3.24|) 
will concentrate in time on the line m y x = m x y. 

As outlined at the end of Section [31 the interesting question to answer is to 
describe the equilibrium profile in the v- variable, to understand if the action of the 
interaction (|4.36[) could produce or not Pareto tails. Since the new variables (v,w) 



4. Fokker-Planck models 



(4.36) 



v* = v + [X(a — P) + p, — fi]w 
w* = (1 — A + ft + (J,) w. 



\w*\ = (l-X)\w\. 



(4.37) 



((1 - A + A + M) 2 ) - (1 - A) 2 + ((/i + m) 2 ) < 1. 
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seem more suitable to describe the problem, let us set g(v,w,t) = f(x,y,t). The 
condition \w\ < v implies that the new density vanishes outside the allowed set of 
values, so that 

(4.38) g(v,w,t) = 0, if \w\ > v. 

On the set \w\ < v, g satisfies the linear Boltzmann equation 



(4.39) 



— / ip(v,w)g(t)dvdw 



(ip(v*,w*) — p(v,w))g(t)dvdw^ 



In the rest of the Section, let us assume that the random variables fi and fl are 
identically distributed, with zero mean and variance a 2 . Thanks to conditions 
(|2.12p . the random variables possess in addition moments bounded of any order. 
Let us expand the smooth function ip(v*,w*) in Taylor series up to order two. We 
obtain 



(ip(v*,w 
A 



ip(v,w))) 



i m dip dtp 

(a - p)w- w— 

dv dw 



2 2 9 V 



1 



{a 



«\2 



,dhp_ 

dw 2 



dv 2 2 
2(q - fi)w 



dw 2 



2 jpp_ 
dvdw 



R(v, w). 



Clearly, R(v, w) denotes the remainder of the Taylor expansion. 

Let us set r = et, and let us consider the situation in which A 
/i — > y/e/J, and ji \/efi. Then h t (v, w, r) = g(v, w, t) satisfies 



eA, while 



— / (p(v 1 w)h e (v,w,T)dvdw 



• i i^ 2 p 1 ., ., 



dip 
dv 



dip 

'7T 
aw 



4* 



i m2 2 9 V 
(a-P) W -—r 

dv 1 



1 

2 

dw 2 



dv 2 
- 2{a-f3)w 



d 2 p 

dw 2 



2 d 2 cp 



dvdw 



h t (r)dvdw+ 



h € {r) + R e (v,w). 



We remark that, by construction, the remainder R e (v,w) depends in a multi- 
plicative way on higher moments of the random variables \J~i[i and v^Ai so that 
R e (v, w)/e <^ 1 for e <C 1 (cfr. the discussion in (3TJH7], where similar computations 
have been done explicitly). 

As e — > h e (v,w,T) — ¥ h(v,w,r) satisfying 



(4.40) 



— / p(v,w)h(v,w,T)dvdw = 
dr J 

w o\ 9 <P \ dL P , 1 2 2°~ 2 <P 
A(a _^_„ Aw _ + _ (TlW _ 



1 2 2 d 2 tp 



2 1 dw 2 

Equation (I4.40j) is the weak form of the Fokker-Planck equation 



(4.41) 



dh 



1 2 2 d 2 h 



2" 1 " dv 2 ' 2"' 2 dw 2 
The coefficients in equation (|4.41[) are given by of = ((fi— fi) 2 ) 
Clearly a\ = a\ = 2a 2 if the errors are uncorrelated. 

The formal derivation of the Fokker-Planck equation (I4.41[) can be made rigorous 
repeating the analogous computations of [23 HZ] , which refer to one-dimensional 
models. The meaning of this derivation is that we allow for small changes in a 



1 2 d 2 (w 2 h) 



dh 

- X(a - 0)w— + A 
dv 



h(r)dvdw. 

d(wh) 
dw 



and of = ((/i+/i) 2 ) 



14 



GIUSEPPE TOSCANI, CARLO BRUGNA AND STEFANO DEMICHELIS 



single trade of goods. Then, in order that a macroscopic change would be visible, 
the system needs to wait a sufficiently long time. We call this procedure quasi- 
invariant trade limit. It is interesting to remark that the balance a 2 /A = C is the 
right one which maintains in the limit equation both the effects of the exchange of 
goods in terms of the intensity A, as well as the effects of the randomness (through 
the variance a). As explained in |47j for the case of opinion formation, different 
balances give in the limit purely diffusive equations, of purely drift equations. 

The Fokker-Planck equation (j4.41j) is reminiscent of all parameters which con- 
tributed to the exchange mechanism. In particular, it contains the values a and /3, 
with a + /3 = I, which are the exponent in the Cobb-Douglas utility function. It is 
remarkable that these parameters appear only in terms of their difference, so that 
the case a = (3 leads to a simplification of the equation. Also, the values of A and 
cr 2 are linked to the differential terms of the equation, the former to the drift, and 
the latter to the diffusion. 

Maybe the most important characteristic of the solution to the Fokker-Planck 
equation, is that it is immediate to recognize that it develops fat tails. This property 
follows by evaluating moments relative to the w variable. Let us set tp(v, w) = 
liul 1 ^', with r > 0. Then, assuming that the solution is rapidly vanishing at 
infinity, one obtains 

— — f \w\ 1+r h(v,w,T)dvdw — 
dr J 

-r<T2 — A^j (I + r) J \w\ 1+r h(v,w,T)dvdw. 

The sign of the right-hand side determines the large-time behavior of the solution. 
If r > 2A/cr| the moment m\+ r = j \w\ 1+r h(v,w,t)dvdw blows up as time goes to 
infinity. On the contrary, in agreement with the behavior of the Fourier transformed 
Boltzmann equation studied in Section |3l if r < 2A/<r| the solution concentrates 
along the line w = 0, that is m y x = m x y. Moreover, since h(v,w,t) = whenever 
w > v, 

(4.42) J v 1+r h(v,w 7 T)dvdw > J \w\ 1+r h(v,w 7 T)dvdw. 

Hence, if r > 2A/cr|, also the principal moment of order r with respect to the v 
variable blows up when time goes to infinity, and the stationary solution, which is 
concentrated along the line w = has Pareto tails. 



5. Final remarks 

In this paper we introduced kinetic equations for the evolution of the probability 
distribution of two goods among a huge population of agents. The leading idea 
was to describe the trading of these goods by means of some fundamental rules in 
prize theory, in particular by using Cobb-Douglas utility functions for the binary 
exchange, and the Edgeworth box for the description of the common exchange area 
in which utility is increasing for both agents. Also, to take into account the intrinsic 
risks of the market, we introduced randomness in the trade, without affecting the 
microscopic conservations of the trade, that is the conservation of the total number 
of each good in the binary exchange. 

The linear Boltzmann equation we introduced in Section |3] and its asymptotic 
Fokker-Planck limit studied in Section SI showed that, despite the microscopic 
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constraints, sufficiently high moments of the solution blow up with time. This is 
in contrast with the behavior of the simplified one-dimensional models for wealth 
distribution expressed by equation (ll.6[) . consequent to trades of type (jl.5l) . where 
the microscopic conservation of wealth in the binary exchange has been shown 
to imply an exponential decay of the stationary profile at infinity [241 139] . thus 
excluding the possibility of Pareto tails. 

The model we introduced in this paper can be generalized in many ways. First, 
the binary trade (|2.11j) suggested by the Edgeworth box and considered in Section 
[5] can be modified in many ways, which continue to satisfy all constraints imposed 
by price theory. One of the possible modifications is given by 

P* = p(l + fj) + A/3(g - p) 
l ''' 4 ' ,) q* =q(l + fi)+\a(p-q). 

Analogously to (12. lip , in (|5.43j) < A < 1, and [i and ft, are random variables with 
zero mean and finite variance, which satisfy conditions (|2.12[) . At difference with 
the previous interactions, here the randomness modifies the point in the Edgeworth 
box proportionally to the percentages of agents (not to their difference). The main 
consequence of this choice relies in the fact that now p = q does not imply p* = q* , 
but only (p*) = (q*) (equality in the mean). This new way of looking at randomness 
reflects heavily on the steady state of the linear Boltzmann equation (|3.24[) . which 
does not concentrate along a well-defined line. 

Second, the strategy of using utility functions and Hedgeworth box for the trad- 
ing is not restricted to an agent-based market in which agents exchange two goods 
only. One can easily generalize the analysis of the previous Sections to model mar- 
kets in which there are N > 2 types of goods to be exchanged, according to (|2.15|) 
for each binary trade. 

Also, one can improve the trading by allowing people to possess different utility 
functions. This can be easily done for example by introducing a Cobb-Douglas 
utility function (|2.10[) in which the exponents a and /3 are positive random variables, 
subject to the constraint 

a + p = l. 

This generalization will introduce a change in the Boltzmann equation, while it 
will induce the same Fokker-Planck limit, with the coefficient a — /3 substituted by 

(a- 13). 
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